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1.  INTRODUCTION 

Turbulence  is  generally  conceived  as  a  collection  of 
eddies  of  many  different  sizes  (Hinze,  1975;  Batchelor, 
1953).  The  "qua si- wavelet”  (QW)  model  discussed  in 
this  paper  is  an  attempt  to  develop  a  mathematical 
representation  for  the  turbulence  that  more  closely  re¬ 
sembles  this  physical  picture  than  Fourier  modes  or  cus¬ 
tomary  wavelets.  Like  customary  wavelets(Farge,  1992; 
Meneveau,  1994),  the  QW  representation  is  based  on 
self-similar  localized  functions.  However,  the  orienta¬ 
tions  and  positions  of  the  quasi-wavelets  are  random, 
and  the  QW  basis  functions  are  not  required  to  be  or¬ 
thonormal  or  to  form  a  mathematically  complete  set. 
Some  other  important  features  of  quasi- wavelets  are: 

•  They  naturally  have  ensemble  statistics  close  to  that 
of  real  turbulence  as  a  consequence  of  the  realistic  basis 
functions. 

•  They  can  simultaneously  provide  Information  about 
scales  of  motion  and  spatial  intermittency.  ‘ 

V  V They  potentially  allow  simplified  models  of  anisotropy 
and  inhomogeneity. 

•  They  can  readily  be  used  to  generate  synthetic  tur¬ 
bulence  fields. 

With  regard  to  this  last  point,  quasi-wavelets  can 
serve  as  a  substitute  for  random  Fourier  modes,  which 
have  previously  been  used  in  many  applications  such  as 
structural  wind  loading  and  simulation  of  wave  scatter¬ 
ing  {Mann,  1998;  Gilbert  et  a!.,  1990).  In  these  applica¬ 
tions,  the  spatially  localized  nature  of  the  quasi  wavelets 
can  be  advantageous  (deWolfe,  1983;  Goedecke  and 
Auvermann,  1997). 

This  paper  is  organized  as  follows.  In  Section  2, 
we  provide  an  overview  of  the  formulation  of  the  QW 
model.  A  relationship  between  the  quasi-wavelet  basis 
function  and  the  energy  spectrum  of  the  turbulent  ve¬ 
locity  fluctuations  is  derived  in  Section  3.  Several  pos¬ 
sible  QW  bases  and  their  corresponding  energy  spectra 
are  discussed  in  Section  4.  In  particular,  a  QW  basis 
function  that  exactly  yields  the  von  Kirmin  spectrum 
is  found.  We  also  consider  possible  QW  models  cor¬ 
responding  to  the  modified  von  K3rm£n  spectrum  of 
Kristensen  et  al.  (1989),  which  includes  the  empirical 
Kansas  spectrum  developed  by  Kalmal  et  al  (1972)  as  a 
special  case.  In  Section  5,  we  present  example  results. 


2.  OVERVIEW  OF  THE  QW  MODEL 

The  main  goal  of  the  QW  model  is  to  represent 
turbulent  fluctuations  by  the  simplest  possible  set  of 
localized  structures  that  resemble  actual  eddies.  In  the 
original  formulation  of  the  model  (Goedecke  and  Auver- 
mann,  1997),  these  structures  were  called  "turbules.” 
Here  we  use  "quasi-wavelet  ”  due  to  the  localized  and 
self-similar  nature  of  the  basis  functions.  For  isotropic 
turbulence,  the  simplest  quasi-wavelet  for  solenoidal  ve¬ 
locity  fluctuations  was  found  by  Goedecke  and  Auver- 
mann  to  be  a  rotating  spherically  symmetric  structure 
given  by 

V  =  V  x  A,  A  =  a2n/  (1) 

where  A  Is  a  vector  potential;  o  Is  the  “size"  or  length 
scale  of  the  quasi-wavelet,  b  is  the  location  of  its  center, 
and  Q  is  Its  angular  velocity  parameter.  4n  an  isotropic 
model,  the  angular  velocity  has  a  uniform,  random  dis¬ 
tribution  over  the  4w  solid  angle.  The  scalar  function 
/,  called  the  QW  envelope  function,  is  any  dimension¬ 
less  localized  function  of  its  argument  |r  —  b|/a.  The 
model  turbulent  velocity  field  results  from  superposition 
of  very  many  such  quasi-wavelets.  For  homogeneous 
turbulence,  each  quasi-wavelet  has  uniformly  random 
center  location  inside  a  chosen  turbulent  volume  F. 

Note  that  the  Individual  quasi-wavelets  do  not  nec¬ 
essarily  satisfy  the  fluid  equations.  We  merely  require 
that  an  appropriate  superposition  of  quasi-wavelets 
must  yield  the  important  statistical  properties  of  the 
turbulence.  In  particular,  we  require  in  this  paper  that 
the  QW  model  yield  correct  or  otherwise  physically  rea¬ 
sonable  spatial  spectra  for  all  ranges  of  the  turbulent 
wavenumber  k. 

Many  different  sizes  an  are  used  in  the  QW  model, 
ranging  from  a\,  the  largest  size  chosen,  to  ajy,  the 
smallest.  Clearly,  a\  corresponds  to  an  outer  scale 
(length  scale  near  the  transition  between  the  energy  and 
inertial  subranges),  and  aN  to  an  inner  scale  (length 
scale  near  the  transition  between  the  Inertial  and  dissi¬ 
pation  subranges).  A  fractal  scaling  of  sizes  is  chosen, 
such  that  0,2/ ai  =  ££3/02  —  . . .  —  /aN_i  =  const. 
The  magnitude  0  of  the  angular  velocity  O  scales  with 
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size  a  such  that  the  characteristic  speed  v  =  Oa  = 
(const.)  (a1/3),  which  is  the  scaling  derived  in  the  Kol¬ 
mogorov  energy  cascade  model  (Batchelor,  1953).  As  a 
result  of  these  chosen  scaling  properties,  the  model  pre¬ 
dicts  the  Kolmogorov  spectrum  in  a  well-defined  inertial 
subrange,  for  any  physically  reasonable  choice  of  QW 
function  /,  if  and  only  If  the  number  density  n  of  the 
quasi-wavelets  of  size  a  scales  like  n  =  (const.) (a-3) , 
i.e.,  the  QW  packing  fraction  na 3  must  be  size  invari¬ 
ant. 

Other  results  in  Goedecke  and  Auvermann  (1997) 
include  physically  reasonable  behavior  of  the  predicted 
spectra  in  both  the  energy  and  viscous  subranges.  In 
fact,  the  QW  model  always  yields  energy  spectra  that 
go  like  k4  for  small  k,  then  transition  to  fc-5/3  for  in¬ 
termediate  k  (the  inertial  subrange),  and  then  fall  off 
much  faster  than  k~^z  for  large  k  (the  viscous  sub¬ 
range  and  beyond).  This  small  k  behavior  is  due  to  the 
fact  that  there  is  a  maximum  QW  size  a\ ,  while  the 
drop-off  at  large  k  is  due  to  the  presence  of  a  minimum 
QW  size  aN. 

3.  QUASI-WAVELET  THEORY  OF  VELOCITY 
SPECTRA 

In  Goedecke  and  Auvermann  (1997),  an  expression 
for  the  velocity,  spectral  tensor  #^(k)  was  derived  from 
a  superposition  of  velocity  quasi-wavelets.  As  shown  in 
Goedecke  et  al.  (2002),  this  expression  leads  to  the  fol¬ 
lowing  relationship  between  the  energy  spectrum  E(k ) 
and  the  scale-invariant,  dimensionless  QW  spectral  func¬ 
tion  F(y ): 


Multiplying  both  sides  of  (2)  by  fe5/3  and  differentiating 
with  respect  to  k  yields 

j-k  [fc^-ECA;)]  =  ala^ika^F^ka,).  (5) 


4.  MODEL  EDDIES 

As  mentioned  earlier,  any  localized  QW  envelope 
/(f)  can  be  normalized  to  yield  exactly  the  Kolmogorov 
(£r5/3)  energy  spectrum  in  an  inertial  subrange.  Only 
the  boundaries  of  the  inertial  subrange  and  the  behavior 
of  the  spectra  outside  it  are  sensitive  to  the  functional 
form  of  /(f).  In  this  section,  we  investigate  several 
possible  QW  envelopes. 

4.1  Von  Karman  eddy 

The  von  Kirm^n  spectrum  is  commonly  used  in  tur¬ 
bulence  modeling.  Application  of  this  spectrum  to  at¬ 
mospheric  turbulence  along  with  appropriate  parameter 
values  is  discussed  by  Ostashev  and  Wilson  (2000).  The 
equation  for  the  energy  spectrum  is 


EvK(k)  — 


CVKalk4Ll 

(l+fe2£2)17/6 


^  55r  (5/6) 


Using  (5),  we  now  find 


FVk{v)  = 
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E(k)  =  <r^a1(fea1)“5/3  /  dy  y14/3  (F(y)f  .  (2) 
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Here,  k  is  the  wave  vector,  is  the  variance  of  one 
of  the  velocity  components,  and  y  =  kat  so  that  the 
integral  extends  over  all  eddy  sizes  a .  The  lower  limit 
on  the  integral  is  actually  kaNt  where  aN  is  the  smallest 
scale  length  used.  We  put  a#  — *  0  here;  this  can  only 
influence  the  behavior  of  E(k )  in  the  viscous  subrange. 
The  QW  envelope  function  /(f)  follows  as  the  three- 
dimensional,  inverse  Fourier  transform  of  F(y): 

m  =  p If  jdZ&yiny),  (3) 

where  £  =  (r-b ) /a  and  f  =  |£|.  Also,  (r  —  b)  is 
the  vector  from  the  uniformly  random  position  b  of  the 
center  of  an  eddy  to  the  general  position  vector  r,  so 
that  f  is  the  scaled  distance  (in  units  of  a)  from  the 
eddy  center.  Because  of  the  spherical  symmetry,  we 
have  in  general  from  (3) 

/(£)  =  ^  f*y  yaia(v0Fiv)-  (4) 


The  von  KSrmin  QW  envelope  function  fvjc( f)  is  de¬ 
termined  by  substituting  (7)  into  (4).  The  integral  can 
be  found  in  tables,  resulting  in 


fvK(§  = 


T(23/6)  0ll1/2  2-11/12 

T(l/3)  Lv\  rV 4r(23/12) 


where  Kv  is  the  modified  Bessel  function  of  the  second 
kind.  The  value  of  a*  is  flexible,  although  since  aj  and 
Lv  are  both  outer  scales,  we  must  have  Lv/ax  ~  1. 
Note  that  in  the  inertial  subrange,  E(k )  =  (55a/l 8) 
62/3  £-5/3^  w^ere  e  is  the  rate  of  dissipation  of  the  tur¬ 
bulent  kinetic  energy  per  unit  mass  and  a  =  0.52  is  an 
empirical  constant.  Matching  this  to  the  von  Kirmin 
energy  spectrum  in  the  inertial  range  kLv  »  1  yields 
a  relation  among  LinaV)  and  e: 


Lv  = 


‘  2T(5/6) 
ay^T(l/3) 
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(9) 


4,2  Generalized  von  Karman  and  Kansas  spectra 

- Kristensen  et  3i^(1989)  proposed  the  following  g€n- 

eral  equation  for  the  one-dimensional  longitudinal  spec¬ 
trum: 

X||()  w(l  +  (kLvM^’  (0) 

with  b(fi)  =  ^r(5/6/i) /T(l/2fi)T(l /3/x) .  The  para¬ 
meter  fj.  controls  the  sharpness  of  the  transition  between 
the  energy  and  inertial  subranges;  ft  =  1  corresponds 
to  the  von  K£rm£n  spectrum.  The  longitudinal  inte¬ 
gral  length  scale  for  this  model  is  £  =  b(fi)Lv.  The 
corresponding  energy  spectrum  is 


Egvk(^) 


Bn 

[(2  —  2/i)  +  (ll/B)(kLv)2^} 


If  /i  =  1,  this  reduces  to  the  von  K5rm£n  expression. 
For  small  k,  the  energy  spectrum  E(k )  for  isotropic 
turbulence  is  expected  to  go  like  k 4  (Hinze,  1975).  We 
note  that  EGvk  goes  like  k4  for  small  k  only  for  /*  =  1 
or  2.  The  empirical  spectra  developed  by  Kaimal  et 
ai  (1972)  on  the  basis  of  the  1968  Kansas  experiment 
correspond  to  fi  =  1/2 ,  for  which  EGVK  is  proportional 
to  k  for  small  k.  For  the  Kansas  spectra  with  /r=  1/2, 
(5)  and  (11)  yield 

*r*/2f.*  _  /40  (LM  (1  +  UL^/Za^* 

27  (i  +  Wai)!4/6 

(12) 

as  the  Fourier  transform  of  the  QW  function  /ka(0- 
This  FKa  diverges  for  y  0,  a  direct  result  of  the 
unusual  behavior  EKa  oc  k  for  small  k .  We  have 
not  attempted  to  find  a  closed-form  expression  for  the 
“Kansas  eddy”  /kaH)*  mainly  because  the  unphysi¬ 
cal  behavior  of  FKa  for  small  y  casts  doubts  that  the 
Kansas  spectrum  can  be  modeled  by  quasi-wavelets. 
Perhaps  this  difficulty  is  to  be  expected,  because  em¬ 
pirical  spectra  can  contain  anisotropic  features  that  are 
not  present  in  the  isotropic  QW  model. 

4.3  Gaussian  and  Exponential  Eddies 

Although  the  von  Kirmin  spectrum  is  simple,  its 
QW  envelope  function  fvK{i)  is  not.  We  now  show 
that  simpler  QW  functions  can  produce  spectra  very 
close  to  the  von  K3rm3n  spectrum.  First,  consider  the 
following  Gaussian-type  envelope,  which  was  used  in 
Goedecke  and  Auvermann  (1997): 

Fa(y)  =  Fg{ 0)e-*2/4  ^  /g(|)  =  7T-^FG(o)e-(\ 

(13) 

where  the  second  equation  results  from  (3).  Substi¬ 
tuting  (13)  into  (2)  and  requiring  that  the  resulting 
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Figure  1:  Comparison  of  various  QW  envelope  functions. 

spectrum  agree  with  the  von  K3rm3n  spectrum  in  the 
inertial  range,  we  obtain 


Fg( 0)  = 


CyK 

2n/6r(i7/6) 


EG{k)  =  -CvK°]  ^(17/6,'fc2af/2) 

Gl  ’  *5/3£j/3  r(17/6)  .  ’  U5J 

where  7 (pfz)  is  the  incomplete  gamma  function. 

We  may  assign  the  ratio  in  Lv/at  in  several  rea¬ 
sonable  ways.  For  example,  we  could  require  the  vari¬ 
ances  of  the  QW  and  von  K3rm3n  spectra  to  be  the 
same.  After  a  lengthy  analysis,  this  leads  to  the  con¬ 
dition  Lv/ai  =0.8379.  Alternatively,  we  could  require 
the  small  k  behavior  of  the  Gaussian  QW  and  von  K3r- 
m 4n  spectra  to  be  the  same.  This  choice  is  discussed 
in  more  detail  in  Goedecke  et  al.  (2002). 

A  second  example  of  a  simple  QW  envelope  con¬ 
sists  of  eddies  with  an  exponential  Fourier  transform 
(exponential  FT  eddies),  given  by 


FB(y)  =  FE(0)e-v^fB($=- 


££.•  <*> 


(i+O 


Substituting  (16)  into  (2)  and  requiring  that  (6)  and 
(2)  agree  in  the  inertial  range,  we  obtain  FE( 0)  and 
the  energy  spectrum  EE(k)  for  these  eddies: 


Fe{  0)  = 


2l7/3CVK 

r(i7/3) 


EE(k) 


CvKoj  1  [7(17/3,2^) 
L  r(17/3) 


Setting  a%  =  we  find  Lvja\  =  0.7238. 


Figure  2:  Comparison  of  energy  spectra  resulting  from  var¬ 
ious  QW  envelope  functions. 


5.  RESULTS 

Various  QW  envelopes  /(f)  derived  in  the  preced¬ 
ing  section  are  shown  in  Figure  1.  The  Gaussian  and 
exponential  envelopes  are  similar,  both  being  flat  at 
the  origin.  The  von  K3rm3n  envelope,  in  contrast,  has 
non-zero  slope  at  the  origin.  Two  versions  of  the  von 
K3rm3n  envelope  are  shown,  one  for  variance  equal¬ 
ing  the  Gaussian  spectrum,  and  the  other  for  variance 
equaling  the  exponential  spectrum. 

In  Figure  2,  we  plot  the  normalized  energy  spectra 
for  the  von  K3rmSn,  Gaussian,  and  exponential  QW 
envelopes.  The  spectra  have  been  normalized  to  have 
the  same  variance.  While  they  are  all  identical  in  the 
inertial  subrange,  in  the  energy  subrange  they  share  only 
the  same  slope  (&4).  The  position  of  the  asymptote 
differs  because  the  integral  length  scale  depends  on  the 
QW  envelope. 

6.  SUMMARY  AND  DISCUSSION 

The  von  K3rm3n  energy  spectrum  of  turbulent  ve¬ 
locity  fluctuations  has  been  widely  used  in  studies  of 
turbulence  and  wave  (acoustic  and  electromagnetic) 
propagation  in  random  media.  In  this  paper,  we  found  a 
QW  envelope  function  that  yields  exactly  the  von  K2r- 
min  velocity  spectrum.  We  also  showed  that  the  QW 
model  has  flexibility  extending  beyond  the  von  K3rm3n 
spectral  model.  In  particular,  it  allows  velocity  spectra 
that  reduce  to  the  Kolmogorov  spectrum  in  the  inertial 
subrange  but  are  adjustable  in  the  energy  subrange. 
This  is  important  because  models  based  on  the  von 
Kirm in  spectrum  sometimes  do  not  agree  well  with  ex¬ 
perimentally  determined  one-dimensional  spectra  in  the 
energy  subrange.  An  objective  in  further  development 
of  the  QW  model  is  to  determine  a  QW  function  and 
corresponding  spectrum  that  yields  the  best  match  to 


experimental  data  in  both  the  energy  and  inertial  sub¬ 
ranges.  This  will  require  a  QW  model  of  anisotropic 
tu rbu fence,  which  we  are  now  attempting  to  formulate" 
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